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Abstract. Verification of the Visible Infrared Imager Radiometer Suite (VIIRS) End-to-End 
(E2E) sensor calibration is highly recommended before launch, to identify any anomalies and 
to improve our understanding of the sensor on-orbit calibration performance. E2E testing of the 
Reflective Solar Bands (RSB) calibration cycle was performed pre-launch for the VIIRS Fight 
1 (FI) sensor at the Ball Aerospace facility in Boulder CO in March 2010. VIIRS reflective 
band calibration cycle is very similar to heritage sensor MODIS in that solar illumination, via 
a diffuser, is used to correct for temporal variations in the instrument responsivity. Monochro- 
matic light from the NIST T-SIRCUS was used to illuminate both the Earth View (EV), via 
an integrating sphere, and the Solar Diffuser (SD) view, through a collimator. The collimator 
illumination was cycled through a series of angles intended to simulate the range of possible 
angles for which solar radiation will be incident on the solar attenuation screen on-orbit. Ideally, 
the measured instrument responsivity (defined here as the ratio of the detector response to the 
at-sensor radiance) should be the same whether the EV or SD view is illuminated. The ratio of 
the measured responsivities was determined at each collimator angle and wavelength. In addi- 
tion, the Solar Diffuser Stability Monitor (SDSM), a ratioing radiometer designed to track the 
temporal variation in the SD BRF by direct comparison to solar radiation, was illuminated by 
the collimator. The measured SDSM ratio was compared to the predicted ratio. An uncertainty 
analysis was also performed on both the SD and SDSM calibrations. 
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1 INTRODUCTION 

VIIRS FI sensor is a cross-track scanning radiometer developed for the NPP (NPOESS Prepara- 
tory Project) mission and is capable of making continuous global observations of the Earth in- 
tended for use in scientific studies of land, ocean, and atmosphere [1-3]. The design of VIIRS 
was heavily influenced by heritage sensor MODIS [4]. VIIRS RSB consists of 14 spectral bands 
covering a range of ~ 0.4 — 2.3 firn with resolutions at nadir of 375 m (3 bands) and 750 m 
(11 bands). In order to maintain the reflective band calibration on-orbit, VIIRS was equipped 
with three on-board calibration devices, the Solar Diffuser (SD), the Solar Attenuation Screen 
(SAS), and the Solar Diffuser Stability Monitor (SDSM). 

VIIRS FI sensor will rely on regular solar observations on-orbit for RSB calibration (similar 
to the heritage sensor MODIS [5]). These observations can occur once per orbit over the South 
pole for approximately three minutes. During these observations, solar illumination, scattered 
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off the SD, will be used to correct for temporal variations in the instrument radiometric respon- 
sivity over the sensor lifetime. However, the SD is fabricated from space-grade Spectralon, the 
performance of which is known to degrade over time from exposure to the sun. The degradation 
of the SD BRF will be tracked by the SDSM, a ratioing radiometer which relates near simulta- 
neous measurements of direct solar illumination to solar illumination diffusely reflected off the 
SD. It is important that the RSB on-orbit calibration cycle be performed during pre-launch test- 
ing, in order to identify any issues and/or anomalies as well as test the validity of the calibration 
cycle under controlled conditions. 

E2E testing of the VIIRS FI RSB calibration cycle was performed in March 2010 at the 
Ball Aerospace facility in Boulder CO. VIIRS had been successfully integrated into the NPP 
spacecraft two months earlier. The NIST T-SIRCUS (Traveling Spectral Irradiance and Ra- 
diance Responsivity Calibrations using Uniform Sources) [6] was transported to Boulder to 
act as a source for E2E testing (as well as spectral testing [7]). The T-SIRCUS was used to 
illuminate both the EV (through an integrating sphere) and SD view (via a collimator) using 
monochromatic light and was monitored in real-time. In addition, both paths into the SDSM 
were simultaneously illuminated. 

E2E testing was not part of the baseline validation testing of VIIRS FI sensor performed 
at the sensor subcontractor (Raytheon El Segundo). The E2E test was conducted during NPP 
spacecraft level testing to investigate whether this type of test setup could validate the VIIRS 
RSB calibration cycle pre-launch. Further, the use of tunable lasers in the calibration is in the 
exploratory phase. As such, this work is focused as much on the application of this technology 
to the calibration of VIIRS as to the results of the analysis. 

This paper will compare the sensor responsivities derived from T-SIRCUS data for both 
views in order to assess the on-orbit SD calibration cycle as well as the SDSM calibration. 
Section 2 will review the design of the on-board RSB calibrators. Section 3 will provide a 
brief overview of the T-SIRCUS and a description of the test setup. Section 4 will establish 
the analysis methodology, and Section 5 will discuss the results of this analysis including the 
measurement uncertainty. 

2 RSB ON-BOARD CALIBRATOR OVERVIEW 

2.1 Solar Diffuser and Solar Attenuation Screen 

VIIRS on-orbit calibration is reflectance based with reference to the SD composed of space- 
grade Spectralon. The SD Bi-directional Reflectance Factor (BRF) was characterized pre- 
launch by the sensor subcontractor (then Raytheon SBRS) [8], The BRF was measured at 
six wavelengths (0.4, 0.5, 0.6, 0.7, 0.9, and 1.7 /jm), nine incident angles designed to replicate 
the range of possible on-orbit solar observations, and two angles of reflectance [corresponding 
to the Rotating Telescope Assembly (RTA) and SDSM views]. 

The SD views the sun through the SAS. Unlike MODIS, VIIRS SAS is fixed; consequently 
only one radiance level is available for on-orbit calibration per band. The vignetting function 
of the SAS was measured pre-launch by the sensor subcontractor [8]; at normal incidence, 
the transmission is ~ 0.13. In addition, no aperture door was incorporated into VIIRS design 
(in contrast, MODIS has an aperture door). As a result, the SD will be exposed to greater 
potential degradation of its BRF on-orbit than MODIS. Improvements in the design and pre- 
launch characterization of VIIRS (particularly with regard to the SDSM) will enable accurate 
tracking of this degradation. 

2.2 Solar Diffuser Stability Monitor 

Degradation in the SD BRF from ~ 0.41 — 0.94 // m will be tracked on-orbit by the SDSM. The 
SDSM is a ratioing radiometer that compares direct solar illumination with light transmitted by 
the SAS and then diffusely reflected off the SD. The direct solar illumination is viewed through 


an attenuation screen, the transmission of which was characterized pre-launch by the sensor 
subcontractor (then Raytheon SBRS) [8]. The SDSM contains an integrating sphere embedded 
with eight detectors filtered to match VIIRS moderate resolution bands in the visible and near 
infrared (with the exception of the eighth detector at 935 nm, which has no corresponding VIIRS 
band). A fold mirror directs light into this integrating sphere from one of three views: solar, SD, 
or dark (the dark view is used as a zero reference for the detectors). A stepping motor rotates 
this fold mirror through these positions every three scans repeatedly during an RSB calibration 
cycle. 

The SDSM collects data during multiple on-orbit RSB calibrations over the sensor lifetime. 
A time-series of the ratios of the solar and SD signals will be used to track the SD BRF degra- 
dation, and to correct the responsivity on-orbit. Due to the fact that the same SDSM optics, 
detectors, and electronics are used for both SDSM solar and SD paths, any slow changes in 
their performance will not effect the quality of the SD BRF degradation monitoring. 

3 TEST OVERVIEW 

3.1 T-SIRCUS 

NIST has developed a laser-based facility known as SIRCUS for use in radiance and irradiance 
responsivity calibrations [6]. A portable version known as T-SIRCUS (or Traveling SIRCUS) 
was transported to Ball Aerospace in Boulder CO during March 2010 for VIIRS spacecraft level 
testing. SIRCUS consists of a series of tunable lasers; for E2E testing, a Ti:sapphire laser was 
used in continuous wave mode to produce monochromatic illumination at 742 and 852 nm (and 
frequency doubling was used to reach 442 nm). 

The T-SIRCUS was used to feed two sources via fiber optics: an integrating sphere and a 
collimator. The integrating sphere was positioned in the EV at a scan angle of ~ —0.2 degrees 
(VIIRS nadir view is defined as 0 degrees). The sphere was approximately 30 inches in diameter 
with an exit aperture of about 10 inches. The sphere was mounted on a scissor jack with its exit 
aperture approximately 4 ft from the telescope aperture at nadir. Two radiance monitors were 
used to track the radiance output from the sphere in real-time. The monitors were positioned 
just off the optical path, one above and one below. 

The collimator was positioned to fully illuminate both the SD through the SAS and the 
solar entrance aperture of the SDSM. The exit aperture of the collimator was about 25 inches in 
diameter. The collimator was also mounted on a scissor jack, by which the angle of the incident 
light on the SD was varied. A total of seven collimator positions were performed covering five 
different angles, listed in Table 1 in the order performed. Note that the first, fourth, and fifth 
positions are repeated measurements of roughly the same angle. An irradiance monitor was 
positioned near the center of the collimator exit aperture, from which real-time tracking of the 
irradiance was performed. 

3.2 Collimator Uniformity 

The collimator beam was known to be spatially nonuniform. Different portions of the collimator 
beam illuminated the irradiance monitor, SD, and SDSM views. As a result, correction factors 
relating the collimator illumination at VIIRS SD view as well as the SDSM SD and solar views 
to the irradaince monitor view were calculated. 

After most collimator positions, the collimator was repositioned to illuminate a large Spec- 
tralon panel at roughly the same angle as was incident on VIIRS. A CCD camera was used to 
photograph the Spectralon panel. From these images, the correction factor was determined. 


Table 1. Collimator positions measured during E2E testing. Angles are relative to VIIRS coor- 
dinate system. 


Position 

Date 

Declination 

Azimuthal 

Wavelengths (nm) 

1 

3/19/2010 

22.52 

16.31 

742, 852 

2 

3/19/2010 

13.64 

16.87 

742, 852 

3 

3/19/2010 

30.09 

16.34 

742, 852 

4 

3/19/2010 

22.37 

16.44 

742, 852 

5 

3/20/2010 

22.38 

16.52 

742, 852 

6 

3/20/2010 

30.44 

15.21 

742, 852 

7 

3/20/2010 

13.63 

15.07 

442, 742, 852 


3.3 Test Description 

For E2E testing, VIIRS was operated in diagnostic mode, during which 2048 unaggregated 
moderate resolution (750 m at nadir) EV samples were recorded covering scan angles from 
about -18.2 to 18.2 degrees. This corresponds to a subset of the full operational mode EV data. 
Each sample has an angular resolution of ~ 0.017776 degrees. The integrating sphere was 
positioned in the EV and fully illuminated samples between scan angles of ~ —1.0 to ~ 0.7 
degrees. During the course of each scan, 48 moderate resolution samples were also recorded in 
each of the three calibration views: Space View (S V), SD, and On-Board Calibrator Blackbody 
(OBC BB). In addition, VIIRS was operated in fixed high gain mode, which restricts all dual 
gain bands to their low radiance settings. 

VIIRS data for E2E testing was recorded in a three data collection cycle. During the first 
collection, the SD and SDSM were illuminated, while the EV was dark. For the second collec- 
tion, the EV was illuminated, while the SD and SDSM views were dark. For the last collection, 
both the EV and SD views were dark. Each of these collections contained 128 VIIRS scans 
(roughly 3.8 minutes). For most collimator positions, this cycle was conducted twice per mea- 
sured wavelength (except position 5 when it was conducted only once). Data collections were 
performed with the SIRCUS output at 442, 742, and 852 nm (442 nm was only measured once 
during collimator position 7); this corresponds to VIIRS bands M2, M6, and M7 (see Table 2). 

In addition, during the first collection in each cycle, the SDSM mirror was rotating and data 
was recorded. The SDSM detectors corresponding to the bands measured are listed in Table 2. 

Table 2. VIIRS bands (and corresponding SDSM detectors) used in E2E testing. 


Band 

Center A (nm) 

Bandwidth (nm) 

SDSM Detector 

M2 

445 

18 

2 

M6 

746 

15 

6 

M7 

865 

39 

7 


4 METHODOLOGY 
4.1 SD Calibration 


In order to test the SD calibration cycle, the calculated responsivities (g) from both the EV and 
SD paths were compared, or 

gsD ,,, 

v = — ■ (!) 

9ev 

Here the detector responsivity is defined as the ratio of the detector response (background sub- 
tracted digital response) to the at-sensor radiance. Ideally, the detector responsivities for the EV 
and SD paths should be equal. As a result, the calculated ratio will depend on how well the 
components in the optical path were characterized as well as knowledge of the input radiance 
(or irradiance). 

The responsivity for the EV path is the detector response divided by the product of the 
integrating sphere radiance and the Response Versus Scan angle (RVS) at the Half Angle Mirror 
(HAM) Angle of Incidence (AOI) of the integrating sphere, or 
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The RVS is the variation in instrument response to HAM AOI and was described in [9]. For this 
test, the EV HAM AOI was ~ 36.3 degrees. The EV radiance {Lev) was determined from the 
average radiance retrieved from the two monitors on the integrating sphere. 

The responsivity for the SD path is the detector response divided by the product of the 
collimator irradiance, the SAS transmission, the BRDF of the SD at the RTA reflectance angle, 
the cosine of the projection angle onto the SD, and the RVS at the SD scan angle, or 
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For VIIRS, the RVS has been normalized to the SD HAM AOI of ~ 60.2 degrees (i.e. II V Ssd = 
1). The SAS transmission is defined as [8] 


tsas(S, (ft = 0.1258 [1 - 0.1538 tan(5)] [1 - 0.04746 tan(<£)] , ( 4 ) 

where <5 and 0 are the collimator declination and azimuthal angles, respectively. The reflectance 
off the SD is assumed to be Lambertian in the neighborhood of the RTA reflectance angle 
(i.e. BRFrta = kBRDFrta )■ The BRF was measured at component level testing for nine 
incident angles, two reflectance angles, and six wavelengths. The BRF was fit to the following 
function of declination and azimuthal angles for each measured wavelength and reflectance 
angle [10]: 

BRF(5, <fi, A) = co (A) + ci(A)<5 + C 2 (A)<j> + C 3 ( A)<5 2 + C 4 (A )</> 2 + C 5 ( A)<5<^>. (5) 


The fit coefficients for the RTA reflectance angle are listed in Table 3. The final BRF is deter- 
mined by interpolating Eq. (5) between wavelengths. The cosine of the projection angle 9 rd is 
the dot product of the sun vector and the SD normal vector [8], or 


cos (9sd) = n sun ■ usd 

where n SD = (0.29724, -0.21860, 0.92944) and 

1 


i/l + tan 2 (S) + tan 2 ((ft 


[1, — tan(<)>),tan(<5)]. 


(6) 

(7) 


As discussed in Section 3.2, the collimator output is known to be spatially nonuniform. As 
a result, a correction factor ( 71 ) relating the collimator footprint on the SD to the irradiance 
monitor position is included, or 

Esd — 7l Emom (8) 

where E mon is the irradiance monitor output. The derivation of the collimator uniformity cor- 
rection factor will be described in Section 4.3. As a result, Eq. (3) becomes 

f = ndnsp 
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Table 3. SD BRF fit coefficients for RTA reflectance angle [10]. 


Wavelength (nm) 

co 

Cl 

C 2 

C 3 

c 4 

C 5 

400 

0.999406 

0.000947 

-0.002517 

-0.000004 

- 0.000010 

0.000026 

500 

1.008981 

0.000454 

-0.002655 

0.000004 

-0.000009 

0.000032 

600 

0.998804 

0.001787 

-0.002700 

-0.000044 

-0.000009 

0.000032 

700 

0.966439 

0.005746 

-1.002742 

-0.000157 

-0.000009 

0.000032 

900 

0.987594 

0.002276 

-0.002797 

-0.000052 

- 0.000011 

0.000025 

1700 

0.962194 

0.005435 

-0.002978 

-0.000153 

- 0.000011 

0.000032 


Table 4. SD BRF fit coefficients for SDSM reflectance angle [10]. 


Wavelength (nm) 

co 

Cl 

C 2 

C 3 

C 4 

C 5 

400 

0.906011 

0.005525 

- 0.001012 

-0.000064 

- 0.000020 

0.000014 

500 

0.976150 

-0.002411 

-0.001104 

0.000161 

- 0.000022 

0.000013 

600 

0.925660 

0.003218 

-0.001126 

-0.000004 

- 0.000022 

0.000014 

700 

0.921604 

0.003537 

-0.001184 

-0.000008 

-0.000023 

0.000013 

900 

0.933213 

0.000525 

-0.001686 

0.000100 

-0.000024 

0.000039 

1700 

0.933563 

0.000067 

-0.001312 

0.000096 

-0.000027 

0.000007 


4.2 SDSM Calibration 

The measured SDSM ratio is defined as follows: 


Rmeasured 
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The Rmeasured is determined every three scans and is trended over a particular calibration cycle 
as well as across calibration cycles. This allows for tracking of the temporal changes in the SD 
BRF. 

The irradiance reaching the SDSM integrating sphere aperture using the SD path is 


Esdsm-sd = 72 E mon TsAS SDSAI cos(0,s£>)7rsin"WO, 


( 11 ) 
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where tsas and cos {Osd) are defined in Eqs. (4) and (6), respectively. The SD BRF is de- 
scribed by Eq. (5) using the fit coefficients for the SDSM reflectance angle listed in Table 4. 
The factor of n sin 2 (ip) corresponds to the solid angle of the entrance cone on the SDSM for 
the SD view, where the half angle of the cone, is 7.78 degrees [8]. 72 E mon is the measured 
irradiance of the collimator (including uniformity correction factor 72). 

The irradiance illuminating the SDSM integrating sphere aperture using the solar path is 

EsDSM-solar = lzE m o n TSDSM ■ (12) 

tsdsm is the transmission factor for the SDSM screen. This transmission includes the cosine of 
the projection angle onto the SDSM screen. In addition, because the detectors are located inside 
the SDSM integrating sphere, each detector has a direct view of the reflection of the image on 
the inside of the sphere. This results in each detector observing a slightly different radiance 
based on its location inside the integrating sphere [8], Note that the collimator uniformity 
correction for the SDSM SD view is different from the correction for the solar view of the 
SDSM. 

The theoretical SDSM ratio is the ratio of Eqs. (11) and (12), or 

^calculated = COs{0 S d) B RFs DS M Sm 2 (/0). (13) 

73 T SDSM 

The source irradiance cancels (excepting the uniformity corrections) and what remains is the 
ratio of the transmission factors in each optical path. Note that the transfer function from the 
SDSM integrating sphere aperture to the SDSM detector is common to both paths and would 
also cancel in the ratio. 

4.3 Collimator Uniformity Correction 

Table 5. Collimator uniformity correction for the SD calibration versus collimator position. 


Wavelength (nm) 

1 

2 

3 

4 

5 

6 

7 

742 

0.967 

0.965 

0.966 

0.964 

0.958 

0.959 

0.959 

852 

0.940 

0.937 

0.933 

0.936 

0.924 

0.926 

0.933 


4.4 Uncertainty Analysis 

The following uncertainty analysis uses the standard error propagation to combine the individual 
uncertainties involved in the present experiment to a total uncertainty [11], 

The measurement error in the EV responsivity is given by the standard error propagation, or 

u 2 {g E y) _ u 2 {dn E y) u 2 {RVS ev ) u 2 {L EV ) 

g E y dn 2 EV RVS 2 ev L'e V 

The SD responsivity uncertainty is also given by the standard error propagation, or 

u 2 {gsp) _ u 2 {dn S p) u 2 {RVS S p) u 2 (7i) u 2 {E S p) 

dSD drigp RVS ED 7i Eg D 
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Since most of these factors were measured independently [tsas, BRFrta, Lev , Esd, and 
cos(Osd)], it is assumed that these factors have no correlations with any other factor. In addi- 
tion, although dnsD, driEV , RVSev , and RVSsd were all measured using VIIRS detectors, 
they were measured separately and therefore are also assumed to be uncorrelated with the other 
error sources. The measurement error in the responsivity ratio is then 


u 2 {rj) _ u 2 (g S p ) u 2 (g E y) 

V 2 Ssd 9ev 


(16) 


The error in the calculated SDSM ratio is 


^ (-^ calculated ) 


td 2 

^ calculated 


u 2 {tsas ) u 2 (tsdsm) M 2 [ cos ( 6 > gp )] 
t sas t sdsm cos 2 (9sd) 
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Since all of these factors in R ca i C uiated were measured independently, it is assumed that there 
are no correlations between the uncertainties of the various factors. 


5 ANALYSIS RESULTS 

5.1 SD Calibration 

5. 1. 1 Data Reduction 

The integrating sphere aperture covered an arc of about 2 degrees of scan angle centered roughly 
on nadir. From this region, 90 moderate resolution (750 m) samples were extracted (samples 
970-1059). These 90 samples were averaged for every scan (ideally there are 128 scans per data 
collection) and for each VIIRS detector (16 detectors per VIIRS moderate resolution band). 
In addition, the 48 samples extracted from both the SD and OBC BB views were similarly 
averaged (note that calibration view data is reported in 14 bits while EV data is automatically 
truncated from 14 to 12 bits by VIIRS; in this analysis all calibration view data is first truncated 
to 12 bits to remove any potential bias). For each average, those points more than three standard 
deviations from the mean were excluded. The OBC BB was maintained at ~ 292 K throughout 
the testing; as a result, the OBC BB data was used as a zero reference for both the EV and 
SD data. The resulting scan dependent driEV and dnsD were substituted into Eqs. (2) and (9), 
respectively. 

The integrating sphere radiance and collimator irradiance were tracked in real-time by mon- 
itors (two monitors on the sphere and one monitor on the collimator). The output of each 
monitor was provided in volts, and was converted to either radiance or irradiance (the average 
of the two radiance monitors is used here) using predetermined monitor calibration coefficients 
provided by NIST. Output from each monitor was obtained about every 10 s. Linear interpola- 
tion was used to acquire the radiance or irradiance at the particular VIIRS time stamp associated 
with each scan. The resulting scan dependent Lev and Bsd were input into Eqs. (2) and (9), 
respectively. 

Using Eqs. (2) and (9), the EV and SD responsivities were determined on a scan by scan 
basis. Since the EV and SD data was recorded at separate times, the scan averaged EV and SD 
responsivities were calculated for each data collection. These scan averaged responsivities were 
then substituted into Eq. (1), from which the responsivity ratio was determined for each three 
data collection cycle. The average of the two collection cycles per collimator position was then 
calculated. 

The following sections will only discuss the results for HAM side A, but both HAM sides 
yield consistent results. 
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Fig. 1. M6 scan dependent responsivities (both EV and SD) from collimator position 7, mea- 
surement 2. 

5.1.2 Data Quality 

A number of data quality checks were performed during the processing of the data for the SD 
calibration. Listed below are the major findings of those checks. 

Some of the data collections did not contain the full expected 128 scans of VIIRS data 
(collections at positions 3 and 5, both at 852 nm). For those collections, the data was retained 
for those scans that appeared to be valid data (i.e. they were consistent with the data collections 
in which all scans were present). 

The signal observed in the SD view when the EV was illuminated via the integrating sphere 
was negligible; in contrast, the response reported in the EV when the SD was illuminated via 
the collimator was between 2 - 3 dn for most cases. This is the possible result of scatter from 
the collimator reaching the EV, and as such has no real effect on the calibration. In addition, the 
data from the dark collections was analyzed and the resulting dn was negligible for all cases. 

The standard deviation of the EV dn was determined over scans at each sample and also 
over samples at each scan (similarly for the SD dn). For the SD view, the two methods were 
comparable. However, for the EV, the standard deviation over samples was lower than the 
standard deviation over scans. This indicates that the integrating sphere had higher spatial 
stability than temporal stability (and that the measured SNR was closer to the true sensor SNR 
using the first method). 

The stability of the sources was also tracked via their respective monitors. The two monitors 
on the integrating sphere show some variation over time (they are both individually stable to 
within about 1.1 % at 742 nm and 2.1 % at 852 nm); in addition, the two monitors vary in 
approximately the same manner. However, the monitors are offset from each other by about 2 
% for 742 nm and 0.2 % for 852 nm. The collimator monitor recorded temporal variation of up 
to about 2.9 % (3.3 %) at 742 nm (852 nm) over a particular data collection. 

The EV dn and monitor radiances have similar temporal trends, and as such the calculated 
EV responsivity at each scan is generally consistent overtime. In contrast, the SD dn and colli- 
mator irradiance have different temporal dependencies, with the result being that the calculated 
SD responsivity is not necessarily constant over scans; an example of this behavior is plotted in 
Figure 1 which shows the scan by scan EV and SD responsivities for M6 collimator position 7 
(second data collection). The inability of the collimator monitor to effectively track the incident 
irradiance is a major source of uncertainty for the SD responsivity. 




5.1.3 SD Uncertainty Analysis 

The error propagation for the SD calibration was defined by Eq. (16). The individual uncer- 
tainties that are used in Eq. (16) are described in this section. In general, each uncertainty 
contributor is the combination of precision (random) and accuracy (biases) errors. 

The error in the SD and EV dn is a combination of the precision and accuracy uncertainties. 
The precision error is the standard deviation of the mean over all samples and scans in a given 
data collection. The accuracy uncertainty accounts for any biases between the SD or EV and 
the OBC views (since any biases common to all sectors cancel in the background subtraction). 
The known biases between sectors are either accounted for in the processing (truncation of cal- 
ibration views to 12 bits before processing begins) or not applicable (biases due to aggregation 
or auto/fixed gain differences); as a result, no accuracy uncertainty is included here. 

The RV S uncertainties were determined during ambient phase testing at the instrument level. 
The band average values used in this work are 0.000570 [M6] and 0.000558 [M7], and are used 
for both SD and EV angles. These values are the Root Sum Square (RSS) of fitting and measure- 
ment uncertainties [9]. The uncertainties associated with the EV radiance are the combination of 
precision and accuracy uncertainties. The precision error is determined by calculating the stan- 
dard deviation of the mean over all radiance monitor values for a particular data collection (both 
monitors). The accuracy error is the average offset between the two radiance monitors over a 
data collection. The total radiance uncertainty is the RSS of the precision and accuracy errors. 
The error in the SD irradiance is the standard deviation of the mean over all irradiance monitor 
values in each data collection (precision error). The accuracy error was undetermined. The er- 
ror in the SAS transmission and the SD BRF were determined by the sensor subcontractor [8], 
and are 0.24 % and 1.09 %, respectively. The cosine of the projection angle is determined from 
alignment measurements made using theodolites and Eq. (6). The measurements are assumed 
to be known to 0.01 degrees. That uncertainty is propagated through Eq. (6). The uncertainty 
in the collimator uniformity correction is currently in progress and is left for future work. The 
final responsivity uncertainties are determined by Eqs. (14) and (15). Then, the uncertainty in 
the responsivity ratio at a particular data collection cycle is determined from Eq. (16). 

5.1.4 SD Calibration Results 

The derived M6 and M7 responsivity ratios for each detector and collimator position are shown 
in Figure 2. In the case of M6, all collimator positions are within ~ 2 % of one. In addition, 
the responsivity ratio tends to increase over detectors for all collimator positions (particulary 
positions 1 and 7). For M7, all positions are within ~ 4 % of unity, with the exception of 
position 7 which is between 4 and 6 % below one. Here the detector variation is again roughly 
increasing, but more uniformly over collimator position than for M6. This detector trend likely 
results from detector dependence in the collimator uniformity correction (the RTA footprint 
on the SD is slightly different for each detector and band, which in turn requires a slightly 
different collimator correction). Although E2E testing was also conducted at 442 nm (M2) for 
collimator postion 7, the digital response was too low to accurately determine the responsivity 
ratio ( dnsD ~ 6 and driEV ~ 15). 

Measurement repeatability was investigated both on short time scales (minutes) and long 
time scales (hours). The collimator was positioned three times at roughly the same angle (po- 
sitions 1, 4, and 5 in Table 1). The interval between positions 1 and 4 was about 6 hours and 
the interval between positions 4 and 5 was about 16 hours. For M6 (see the black, green, and 
yellow curves in the upper plot of Figure 2), the ratios are within 1.0 % of unity for all detectors 
for positions 4 and 5; however, the ratios derived from position 1 are lower (between 0.5 and 2 
% below one). In addition, the detector variation is greater for position 1 . In contrast, M7 shows 
the ratios derived from positions 1 and 5 in good agreement with each other (although between 
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Fig. 2. M6 and M7 responsivity ratios versus VIIRS detector. Collimator positions are given in 
(declination angle, azimuthal angle). 

2 and 4 % below one), while the ratios from position 4 are between 0.5 and 1.5 % below one 
(see the black, green, and yellow curves in the lower plot of Figure 2). 

In terms of short term stability, the E2E test was repeated for each collimator position twice 
within about 30 minutes (except position 5). The data from each individual collection cycle 
(SD, EV, dark) was used to derive a responsivity ratio. The comparison of the short term stabil- 
ity is shown in Figure 3 versus detector. Here the solid (dashed) lines indicate the first (second) 
measurement. For M6, positions 4 and 7 agree well, while positions 1, 2, 3, and 6 show dis- 
crepancies up to 1 %. The ratios deteremined for M7 show agreement to within 0.5 % for all 
positions except postion 2, for which short term repeatability is about 2 %. 

The discrepancies between the repeated measurements are largely understood by examining 
the underlying dn, radiance, and irradiance. In Table 6, the dn, radiance, and irradiance are listed 
for M6 detector 9 at each collimator position and measurement repeat (averaged over scans). 
For the EV, the dn increases when the radiance increases; however, this is not always the case 
for the SD measurements. Collimator positions 1, 2, and 7 actually have lower irradiance for 
the measurement with higher dn. This result also holds when examining the SD on a scan basis. 
Thus the discrepancies in the repeatability measurements appear to be largely related to the 
accuracy of the irradiance monitor output. 

The band averaged uncertainty estimates based on Eq. (16) are listed in Table 7 for all 
collimator positions. The BRF uncertainty is the largest contributor to the uncertainty for M7 
(at 1.09 %). For M6, the largest contributor is the accuracy error for the EV radiance, which 
is around 1 .7 %, follwed by the BRF uncertainty. Because the main sources of error were not 
detector dependent, the varition of the uncertainty with detector is small for both bands. Note 
that the collimator uniformity correction uncertainty was not included; it is expected that this 
uncertainty will be a major contributor. 

Figure 4 shows the band averaged responsivity ratio for each collimator position with the 
uncertainties shown in Table 7. The responsivity ratios are consistent with one for band M6 
for all collimator positions. However, this is not the case for band M7, where only collimator 
position 4 is consistent with one; all other positions are lower (especially position 7 which is 
much lower). 
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Fig. 3. M6 and M7 resposnsivity ratios - short term stability. Collimator positions are given in 
(declination angle, azimuthal angle). 
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Fig. 4. Band averaged responsivity ratios with uncertainties. 


Table 6 . Short term repeatability for M 6 detector 9 (HAM side A), dn, radiance, and irradiance 
are averaged over scans. 


Position 

Measurement 

Esd 

dnsD 

Lev 

dnEV 

V 

1 

1 

4.728 

297.7 

0.1062 

292.7 

0.996 

1 

2 

4.735 

297.4 

0.1061 

292.7 

0.993 

2 

1 

4.701 

286.2 

0.1071 

295.2 

0.986 

2 

2 

4.729 

283.7 

0.1061 

292.7 

0.991 

3 

1 

4.674 

290.2 

0.1053 

290.6 

0.997 

3 

2 

4.612 

287.6 

0.1055 

291.3 

1.003 

4 

1 

4.747 

291.2 

0.1076 

296.7 

1.008 

4 

2 

4.735 

291.4 

0.1069 

294.8 

1.008 

6 

1 

5.924 

354.6 

0.1286 

355.2 

1.005 

6 

2 

6.006 

358.1 

0.1319 

364.1 

0.999 

7 

1 

6.089 

369.1 

0.1286 

354.4 

0.986 

7 

2 

6.084 

369.9 

0.1278 

351.9 

0.986 


Table 7. Responsivity ratio uncertainties in percent. 


Uncertainty 

1 

2 

3 

4 

5 

6 

7 

M 6 1 /( 77 ) 

2.19 

2.10 

2.11 

2.11 

2.08 

2.08 

2.07 

M7 1 /( 77 ) 

1.20 

1.19 

1.20 

1.20 

1.20 

1.21 

1.20 


5.2 SDSM Calibration 

5.2.1 Data Quality and Reduction 

Five SDSM samples were recorded every scan for each detector in the SDSM integrating sphere. 
These five samples were averaged for every scan. The averaged DN from each three scan cycle 
was used to determine the ratio in Eq. (10). This SDSM ratio was trended versus time for 
each data collection and collimator angle. In addition, the average of the ratios for each data 
collection was determined. As the SDSM was illuminated with monochromatic light, only the 
detector corresponding to the input wavelength was analyzed. As with the SD data, some data 
collections were missing scans; again, the data was retained for those scans that appeared to be 
valid data. The collimator uniformity corrections (72 and 73) were not applied in this work. 

5.2.2 SDSM Uncertainty Analysis 

The error propagation for the SDSM calibration was defined in Section 4.4. The uncertainties of 
the measured SDSM ratio as well as the individual uncertainties that are used in the calculated 
SDSM error are described in this section, along with the total uncertainty in the calculated 
SDSM ratio. 

Due to the large variability in the SDSM measurement over any given data collection (as 
shown in Figure 5), the error on the average SDSM ratio per data collection is taken to be the 
standard deviation of all individual SDSM ratios within that data collection. 


The error in the SAS transmission, the SDSM screen transmission, and the SDSM BRF 
were determined by the sensor subcontractor [8], and are 0.24 %, 0.67 %, and 1.09 %, re- 
spectively. The uncertainty in the cosine of the projection angle is determined from alignment 
measurements made using theodolites. The measurements are assumed to be known to 0.01 
degrees. That uncertainty is propagated through Eq. (6). The SDSM entrance cone angle is also 
assumed to be known to 0.01 degrees; the uncertainty is propagated to the solid angle of the 
entrance cone. The determination of the uncertainties for the collimator uniformity correction 
factors is in progress; for the purposes of this work, they are not included in the final results 
presented below. The final calculated SDSM ratio uncertainties are defined by Eq. (17) for each 
data collection. 

5.2.3 Results 

The measured SDSM ratios [as defined by Eq. (10)] are plotted in Figure 5. Each color rep- 
resents a different collimator position (solid and dashed lines indicating first and second data 
collections at a particular collimator position). The ratios are between 1.1 and 1.3 for 742 nm 
and between 1.15 and 1.35 for 852 nm. In addition, the measured ratios fluctuate by up to 5 
% over roughly 3.8 minutes (the variaition in the solar view is generally larger than in the SD 
view). On MODIS, the SDSM also exhibits large variation in any given detector over a typical 
solar observation; fortunately, all the SDSM detectors vary in roughly the same manner. This 
allows the temporal variation to be removed using the response of MODIS SDSM detector 9 
(936 nm). Due to the fact that T-SIRCUS was a narrow band source, only one SDSM detector 
recoded a meaningful signal for a particular data collection; as a result, the temporal trend could 
not be removed in this test using the MODIS methodology. 

The measured ratios, averaged over each data collection, are plotted versus collimator posi- 
tion number in Figure 6 for both 742 nm (upper plot) and 852 nm (lower plot). Black and red 
points indicate the first or second data collections at a particular collimator position (where the 
error bars are ± one standard deviation over that data collection). The blue curve represents the 
calculated ratio at each collimator position. For M6, the calculated and measured ratios agree 
well for all collimator positions. For M7, there is some disagreement between the theoretical 
and measured values as well as some discrepancy between the first and second data collections 
at some collimator positions. 

The measured SDSM uncertainty as defined above was between roughly 0.5 and 2.1 %. The 
calculated error was approximately 1.38% (this number was fairly constant as the major factors 
were not wavelength or angle dependent). The largest contributor was the SD BRF (1.09 %). 
The measurements are consistent with the theoretical values for all positions at 742 nm (with 
the exception of position 3, measurment 1). However, most measurements are not consistent 
with the theoretical results at 852 nm. 

6 CONCLUSIONS 

E2E testing of a limited number of VIIRS FI reflective solar bands was conducted in March 
2010 at Ball Aerospace in Boulder CO using the NIST T-SIRCUS. The E2E test was an explo- 
ration of the use of tunable lasers in the validation of the RSB calibration cycle pre-launch. The 
goal of the testing was to verify that the responsivites determined from the EV and SD paths 
were consistent. In addition, the SDSM measurements were compared to the modeled ratios. 

The EV / SD responsivity ratios were consistent with one for band M6, but inconsistent 
for band M7. Slight upward trending with detector was observed in the responsisvity ratios, 
which likely results from detector dependence of the collimator uniformity correction. The 
SDSM measured and theoretical values are consistent at 742 nm, but inconsistent at 852 nm. In 
addition, the SDSM ratios exhibit some temporal variability. 
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Fig. 5. M6 and M7 SDSM measured ratios - short term stability. 
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Fig. 6. M6 and M7 SDSM ratios (measured and calculated) with uncertainty estimates. 


A large source of uncertainty is in the collimator irradiance, which in turn is heavily depen- 
dent on the collimator uniformity correction. This correction is very complicated and is still in 
progress. It is expected that an updated correction will improve the agreement of the responsiv- 
ity ratio for M7 as well as decrease the detector dependence. In addition, uniformity corrections 
in the SDSM analysis have yet to be determined, which it is hoped will improve agreement for 
the SDSM ratios at 852 nm. 

The E2E test was a preliminary rough validation of the RSB calibration cycle. Improve- 
ments in the test setup should reduce the uncertainties and improve the accuracy of the mea- 
surement (particularly with respect to the collimator irradiance). Future applications of this test 
methodology should expect to accurately validate VIIRS RSB calibration cycle pre-launch. 
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